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Introduction 

Let E be a finite dimensional Euclidean space of dimension d, (d > 2)0- Let 
G ~ SOo(l, d + 1) be the (connected component of) the conformal group of 
E, acting by rational transformations on E. Much interest has been devoted 
to various invariant or covariant objects for this action. There is a natural 
action of G on a space of densities on E. Identifying the densities with 
functions on E, the action is given by 

*x(g)f{x) = K{g-\x)o+ x f{g-\x)) , 

where n(g, x) is the conformal factor (the infinitesimal rate of dilation) of 
the transformation g £ G at x & E, A is a complex parameter and p = i. 
Let Ai, A2, A3 be three complex numbers. A continuous trilinear form C on 



1 The case d = 1 could be treated along the same lines, but there are some differences, 
which would require separate statements. See [15] for a study of this case. 
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C^°(E) x C^°(E) x C^°(E) is said to be conformally invariant with respect 
to ir\ 1 7T\ 2 <8> 7Tx 3 , if, for any three functions /i, /2, /3 G C^°(E) 

£{'K\ i {g)fu'K\ 2 (g)f2,'K\ 3 {g)f2 l ) = C(h,h,h) 

where g is in G (strictly speaking, defined on the union of the supports 
of the three functions). These trilinear forms have been investigated in a 
previous work of the second author in collaboration with B. 0rsted (see [3]). 
Generically, for A = (Ai, A2, A3) in C 3 , there is a unique (up to a multiple) 
such invariant trilinear form. Viewing the trilinear form as a distribution 
on E x E x E, it has a smooth density on the open set 

{(x 1 ,x 2 ,x 3 ) e E x E x E ;xi / x 2 ,x 2 / x 3 ,x 3 / xi} , 

given by 

//3(xi,X2,x 3 ) = \xi - x 2 f 3 \x 2 - a^ 1 !^ - Xif 2 

where j3 = {j3\, f3 2 , P3) is a triplet of complex numbers, uniquely determined 
by A = (Ai, A2, A3) (see (|23]l ). The corresponding distribution Cp is defined 
by meromorphic continuation, and has simple poles along certain planes in 
C 3 . The study of the residues was begun by the second author in [2], and the 
present paper deals with the most singular residues. They are distributions 
supported on the diagonal V = {(x,x,x);x € E}. They turn out to be 
intimately related with covariant bidifferential operators, that is differential 
operators from C^°(E x E) into C^°(E) which satisfy a relation of the form 

D(nx(9)fi ® 7r„(g)f 2 ) = ^{g)D{h ® f 2 ) 

for some (A,/x, v) £ C 3 and for any g E G, fi,f 2 £ C%°(E). Such oper- 
ators have a long history (see [8], [9]), and the most celebrated ones are 
the Rankin Cohen operators, which are holomorphic bidifferential operators 
on the complex upper half plane and covariant under the group PSL 2 (M>). 
In the context of conformal geometry, they were studied by Ovsienko and 
Redou (H7j). See also p], [H]. 

The basic ingredient we use for computing the residues is a Bernstein- 
Sato identity (cf Theorem 14. ip . Recall that for fi, f 2 , f% three arbitrary 
nonnegative polynomials on E, there exists a differential operator B = 
B(xi, x 2 , £3, si, s 2 , S3, d xi , d X2 , d X3 ) on E which is polynomial in the Xj and 
Sj, and a polynomial b on C 3 , such that 
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Such identities exist in general (cf [19], [20] ) , but their explicit determina- 
tion is seldom known. Once such an identity is known, the computation 
of the residues is easy. As a consequence, we find new expressions for the 
covariant bidifferential operators alluded to previously. More general results 
on Bernstein-Sato identities will appear in [JJ. 

The plan of the paper is as follows. Section 1 collects some results on 
distributions supported on a subspace. Section 2 is an elementary approach 
to computing the residues, much in the spirit of Gelfand and Shilov ([TO]). 
In particular it allows to determine the residues along the "first" plane of 
poles. Section 3 exploits the invariance property of the trilinear forms and 
the covariance property of the associated bidifferential operators, making 
connection with the results of [17] , Section 4 is devoted to the Bernstein- 
Sato identity (Theorem I4.ip , which is used in section 5 to give a formula 
for the residues, and as a consequence, a new expression for the covariant 
bidifferential operators (see f)35f) ) . Section 6 presents several remarks and 
perspectives on the subject. 

The second author wishes to dedicate the present paper to Elias Stein, 
on the occasion of his 80th birthday. 

1 Distributions supported on a subspace and trans- 
verse differential operators 

Let E be a finite dimensional real vector space, and let V be a linear subspace 
of E. Let E' be the dual space of E, and let 

V ± = {£ G E' ; £| V = 0} . 

Let u be a distribution on V. The assignment 

C?(E)B<p^(u,tp\ V ) 

defines a distribution on E, the natural extension of u, hereafter denoted by 
u. Clearly Supp(u) = Supp(u) C V. We now characterize the wavefront set 
of u (cf p]). 

Proposition 1.1. Let u be in D'(V), and let u the associated distribution 
on E. Then 

WF(u) = {(ac,0 G Supp{u) x (E'\0); £ G V L or (x,£i v ) G WF(u)} . (2) 
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Proof. Choose a subspace W such that E = V © W. For £ G E', let 
£ = + where £' G VF^ and f" G T/^. Let 99 be in C™(E). Then ^2 is 
a distribution with compact support and its Fourier transform is given by 

Let (xo,£o) be in the set described by the RHS of fl2J). If Co = 0, then 
choose ip such that (it, cp\ v ) 7^ (which is always possible since x belongs 
to Supp(u)), so that F(ipu)(^) cannot decrease rapidly in a conic neigh- 
bourhoud of £ - If 0eo>£o) belongs to WF(u), F{ym){£) cannot decrease 
rapidly on a conic neigbourhood of £q in W 1 - \ 0, a fortiori on a conic neig- 
bourhood of £ hi E' \ 0. Conversely, assume (xo,£' ) does not belong to 
WF(u) and £q 7^ 0. Then, for ip with a sufficently small support near xq 
and for £ in a (small enough) conic neigbourhood of £o> -^"(v 7 ")^) can be 
dominated by Cjv(1 + |£'|) -7V for any integer N. But in a (sufficently small) 
conic neighbourhood of £g one has |£| < C|£'| for some constant C > 0, so 
that T(<pu)(£) is dominated by CW(1 + |£|) _JV . Hence (x ,£o) WF(iI). □ 

To further investigate distributions supported on V, one needs to intro- 
duce normal derivatives. Fix a splitting E 1 = 1/ © VF as above, and choose 
ccordinates wi, W2, ■ ■ ■ , w p on W, which can be regarded as (a partial set of) 
coordinates on E by extending them by on V. Let I = . . . , i p ) be a 
p-tuple of natural integers, let |/| = £1 + £2 + ' • • + V Let Dj be the operator 
(the -D/'s are often referred to as normal derivatives), defined by 



Dup (v) = 




mapping smooth functions on E to smooth functions on V. To any distri- 
bution u on V, one can associate the distribution Dju defined by 

(-lf^Dju,^) = (u,D lV >) . 

Observe that WF(Dju) = WF{u). The inclusion C is obvious, whereas the 
opposite inclusion is obtained by testing against functions (p of the form 

ip(v,w) =x(v)w I ip(w), (3) 

where x £ C^°(V), w 1 = w 1 ^ . . . w V and ip is a function in C^°(W) which is 
identically equal to 1 in a neigbourhood of 0. 
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Now let U be a distribution on E, with Supp(U) C V. The structure 
theorem of L. Schwartz asserts that there exist distributions uj on V such 
that 

I 

where the sum is locally finite. Moreover, the u/'s are unique. 

If all the distributions u\ are given by smooth densities, then from ([2]), 
WF(U) C E x (V 1 \ 0). The converse is true. 

Proposition 1.2. Let U be a distribution supported in V , and assume that 

WF{U) C V x (V L \ 0) . 
Then there exist smooth functions ui on V such that, for any <p G C^°(E) 

(U,(p)= / y^u I (v)D I ip(v)dv . 
Jv j 

Proof. By the previous result, U = Djuj, where uj is some distribution 
on V . The assumption on the wavefront set of U, when tested against 
the functions of the form given by (|3|) implies that, for each <i-tuple /, 
WF(ui) C V x (F ± \0), which implies that WF(ui) = by Proposition O 
As the projection onto the first coordinate of the wavefront set is precisely 
the singular support, each u\ coincides with a smooth function on E. □ 

A transverse differential operator D is a mapping form C^°(E) in C^°(V) 
which is given by 

Dip(v) = ^Jai(u)Z)j^(u) , 
i 

where Dj are the normal derivatives introduced earlier, and the a/'s are 
smooth functions on V . The sum is always assumed to be locally finite. No- 
tice that the ai are well determined, again by testing the operator against 
functions of the form given by ([3]). The previous proposition can be refor- 
mulated as : any distribution U supported on a linear subspace V, such that 
WF{U) C V x (V ± \ 0) can be realized as 

(U,<p)= [ D V (v)dv , (4) 
Jv 

for some transverse differential operator D. Moreover, (once a splitting of 
E as V ® W has been chosen) D is uniquely determined . 

Invariance properties of a singular distribution are reflected in the asso- 
ciated transverse differential operator. Here is a special case, fitted for our 
needs. 
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Proposition 1.3. Let U be in D'(E), supported on V. Assume that U is 
invariant under translations by elements of V . Then 

(U,<p) = [ Dip(v)dv , 
Jv 

where D is a transverse differential operator with constant coefficients. 

Proof. Let v be any element of V \ 0, and let X v be the vector field on E 
which is constant and equal to v at each point of E. Then, the invariance 
property of U amounts to the equalities X V U = for any v £ V. Hence, by 
[TT] Theorem 8.3.1 

WF(U) C {(x, G V x (E 1 \ 0),£(v) = 0} . 
As this is valid for any v £ V, 

WF{U) C {(x, £) e V x {V^ \ 0)}. 

By Proposition 11.21 U is given by a transverse differential operator D, i.e. 

(U,(p) = / Dip(v)dv 
Jv 

where D(p(v) = ai(v)Dnp(v). Now, for any vo G V, X Vo commutes with 
any Dj, such that, by integration by parts, 

= (U,X Vo <p) = - [ Y^X Voai (v)DMv)dv . 
Jv i 

Now fix a cf-tuple /, and check this equality on functions of the form Q. It 
yields X vo ai = for any vq £ V and hence aj is a constant. □ 

Remark. All the results of this section could be formulated for distributions 
supported on a closed submanifold. 

2 An elementary approach to residues at poles of 
the second type 

In this section, we consider the standard Euclidean space E = M. d and denote 
the distance of two points x, y G E by |x — y\. Let (3 be a complex number 
and let 

h( x >y) = \x~ yf ■ 
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For IStp large enough, the kernel lp is locally integrable on E x E, so that it 
defines a distribution on E x E. It can be extended meromorphically (as a 
distribution), having simple poles at the points j3 = —d — 2k, k E N (see e.g. 

[ID]). 

Now let (3 = (/Si, ft, ft) be in C 3 . Set, for x^a^)^ £ ^ 
/ /3 (xi,x 2 ,x 3 ) = lf3 1 {x 2 ,X 3 ,)l(3 2 (x 3 ,X 1 )lp 3 (x 1 ,X 2 ) ■ 
For / in C^°(E x E x E) the integral formula 

£(3(f)= / /(xi,X2,X3)/ l3 (xi,X2,X 3 )(i3;idx2(i3;3 

is well defined for §t(/3j)(j = 1,2,3) large enough and can be meromorphi- 
cally continued in C 3 as a distribution on E x E x E. 

Theorem 2.1. The map (3 i — > Cp can be meromorphically extended to 
C 3 , with simple poles along the four families of planes defined by one of the 
following equations 

Px = -d-2k x , p 2 = -d-2k 2 , p 3 = -d-2k 3 (5) 

Pi + fa + ft = "2d - 2k (6) 

where ki,k2,k 3 ,k £ N. 

The analoguous result for the sphere was obtained in a joint work of the 
second author with B. 0rsted [3], and the proof for the flat case requires 
only minor modifications. In fact the flat case and the spherical case are 
related by a stereographic projection already used in [3]. 

A pole is said to be of first type if it satisfies one of the equations ([5]) , of 
second type if it satisfies one of the equations © . 

For poles of the first type, the residues were studied in [2]. We concen- 
trate on poles of the second type. For each 6 N, let Hk be the plane 
defined by 

H k = {(3 e C 3 ; px + fa + ft = ~2d - 2k} . 

In [3] was already proven that, at a generic point (see precise statement 
below) of such a plane, the residue (viewed as a distribution on E x E x E) 
is supported on the diagonal 

V = {(x,x,x)\x £ E} . 

This will be reproved, in a more elementary and explicit way (see Theorem 

B3D- 
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The distribution £p is invariant by diagonal translations, i.e. by map- 
pings t v ,v G E, where 

t v (x!,X2,X 3 ) = (Xi +V,X 2 + V, X 3 + v) . 

Clearly the residue at some pole (3° will have the same invariance. So, for 
f3 a pole of the second type, we are in the situation of Proposition II. 31 One 
has to choose a normal space to T> in E x E x E, and our choice will be3 

W = {(0,y,z),y eE,zeE} . 

The concept of transverse differential operator introduced in the previous 
section can be reinterpreted in this context. A bidifferential operator D is a 
map from smooth functions on E x E to smooth functions on E of the form 

0\I\+\J\(n 

D<p(v) =J2a u (v) QyiQzj (v,v) , 

where / and J are d-tuples of integers, yj (resp zj) are coordinates on the 
first factor (resp. second factor), associated to a choice of a basis of E and 
the a/j's are smooth functions on E. The sum is assumed to be locally 
finite. When the ajj are constant functions (hence all except for a finite 
number), the operator is said to be with constant coefficients. 

Theorem 2.2. Let (3° G U k for some k G N, but such that none of the 
equations ([5]) is satisfied. Then there exists a unique bidifferential operator 
with constant coefficients D^o such that, for f in C^°(E) and g G C^°(E x E) 

Re8(Cp,0°)(f®g)= [ f{x)D 0S g{x)dx . (7) 

Here / ® g stands for the function 

E x E x E 3 (x,y,z) i — ► f(x)g(x,y) . 

Needless to say, as those functions are dense in C^°(ExExE), ([7]) determines 
completely the residue. 

As already observed, the distribution Cp is invariant by any diagonal 
translation. To take adavantage of this, define for ip G C^°(E x E x E) 



$(y,z)= / <p(v,y + v,z + v) dv . (8) 



E 



2 As symmetry among the three variables is broken, from now on, we use (x, y, z) for 
notation of an element in E x E x E. 
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For ip in C%°(E x E x E), the integral converges and defines a function <J> 
which belongs to C^°(E x E). Moreover, the correspondance ip i — >• $ is 
continuous. Notice for further reference that 

$(0,0) = / ip(x,x,x)dx . (9) 
Je 

Lemma 2.1. Assume that > -d, j = 1,2,3 and + /3 2 + /3 3 ) > 

-2d. T/ien, for any ip G C^°(E x E x E) 

I \yf 3 \zf 2 \y-z\^$(y,z)dydz . 

JExE 

Proof. The conditions on the parameter j3 guarantee the convergence of the 
integrals. The equality of the integrals is obtained through the afhne change 
of variables 

v = xi, y = x 2 -xi, z = x 3 - x\ . 

□ 

Let 

S = {(<t, t) e E x E, \a\ 2 + |r| 2 = 1} 

be the unit sphere in E x E, and denote by dfi the Lebesgue measure on S. 
Recall the integration formula in polar coordinates 

[ ®(y,z)dydz = f f ^(ra,rT)dfj,(a,T) r 2d ~ l dr (10) 
Jexe Jo Jt, 

Lemma 2.2. For mC c °°(S) let 

X p (if>) = J \af 3 \rf 2 \a - rf 1 ^, r) dfi(a, r) . (11) 

i) Assume that > — d for j = 1,2,3. Then the integral (|lip is 
convergent and defines a distribution Ip. 

ii) The map (3 i — ► Xg can 6e extended meromorphically to C 3 , 
simple poles along the family of planes given by the following equations : 

= -d- 2kj, j = 1, 2, 3, kj G N . 

Proof. The three subsets of S 

{(ff,r) GE,a = 0}, {(a, r) G S, r = 0}, {(a, r) G S, a = r} 
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are three disjoint submanifolds of dimension d — 1 (hence of codimension d) 
in X. Recalling that we assumed $l(f3j) > —d, for j = 1, 2, 3, the integrals 

! \af z dn(o-,T), I \rf 2 dfJ,(a,T), I \a -T\ h dn(o,T) 

are convergent and hence the integral Xg is convergent by applying a suit- 
able argument involving a partition of unity. This shows i). Similarly, the 
meromorphic extension and the location of poles (also the fact that the poles 
are simple) are classical and can be easily deduced from [10]. □ 

Let $ be a function in C^°(E x E). For r in R, let ip r be the function on 
£ defined by 

Vv(ff,r) = $(n7,rr), (<x,r)e£. (12) 

Then ?/v belongs to C°°(S) and the map (r, $) i — 5- ?/v is continuous from 
E x C£°(.E x E) to C°°(S). 

Lemma 2.3. Assume i/tcrf > -d,j = 1,2,3 and K(/?i + /3 2 + /%) > 

-2d. T/ien, /or any <p e C™(E x E x E) 



' 00 



Cp{ip)= r M-L+Pi+P*+to Zpijjrdr . (13) 

This is just using the formula for integration in polar coordinates. 

Lemma 2.4. Xe£ 7 oe in C£°(M) and assume that 7 is an even function. 
Then the integral 

Is(l) = / r s ~f(r)dr = - \r\ s j(r)dr 

Jo 2 y_ 00 

is convergent for > —1. The map s 1 — > I s (j) can be extended mero- 
morphically to C with simple poles at s = —1 — 2k, k € N. Moreover, the 
residues at the poles are given by 

*»C«7).-l-») = i^(£)"<0). (14) 



For a proof, see [10] . 

Theorem 2.3. Let k G N and Zef /3° = safe/y i/ie following 

assumptions : 

i) ft + /3° + /3° = -2d - 2fc 
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ii) fj i-d- 2N. 

Let ip be a function in (ExExE) and form successively the functions 
$ defined by ([SD and ip r defined by (|12j) . The function (3 i — )• Cp{<>p) has a 
residue at /3° given by 

* , ^^ = f5* + ij(*U^ ,( * ) - (15) 

Proof. Observe that Ipif} r is well defined (lemma 12. 2p and, as a function 
of r is easily seen to be in C£°(R). Moreover, the distribution Z@ is even, 
whereas ip- r (o~,T) = ip r {—a, — r), hence Xg-i/v is an even function of r. Now 
let 7 = Zpipr- Then (fT3|) can be rewritten as Cp{<-p) = I s (l)- Observe that 
2d - 1 + P$ + /3§ + /S§ = -1 - 2k and eventually apply ((HD to conclude. □ 

The expression obtained for the residue (viewed as a distribution on 
ExExE) shows that it is supported by the diagonal V. In fact, if 92 G C°°(5) 
vanishes on a neighbourhood of 2?, then $ vanishes in a neigbourhood of 
(0,0) in E x E, and hence ip r vanishes identically for |r| small enough, 
so that Ip(ip r ) = for \r\ small enough. Hence the residue (evaluated 
against <p) at (3 vanishes. Now, the structure of distributions supported by 
a submanifold is known from Schwartz's theorem. This requires choosing 
at each point (x, x, x) of T> a complementary subspace to the tangent space 
("normal coordinates") at T>. The choice will be 

M"x,x,x = {(x, y, z),y e E,z e E} . 

Moreover, as we are interested in the trilinear form rather than the distri- 
bution, the space of test functions will be restricted to functions of the form 
(/ ® 9) (x, y, z) = f(x)g(y, z), x e E,(y,z) G E x E. 

With this change of point of view and notation, let us write more explic- 
itly (|15p . We will use the following convention for coordinates on E x E : 
for 1 < i < 2d, let 

Vi = yi, if 1 < i < d, Vi = Zi-d, if d + 1 < i < 2d . 

and similarly on S 

Pi = o~i, if 1 < i < d pi = Ti_d, if d + 1 < i < 2d . 

Theorem 2.4. Let (3° satisfy the same assumptions as in Theorem ] 2. 3[ Let 
f G C c °°(£) and geC^(ExE). Then 

Res(Cp{f®g),(3°) = * + J f{x){D^g){x)dx , (16) 
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where D^o is the bidifferential operator with constant coefficients given by 

d 2k g 

l<h,i 2 ,...,i 2k <2d 11 12 %2k 

where 

a h ,i 2 ,...,i 2k (P°) = J p h . . . Phk \af°\Tf°\a - rf°dn(a,T) . (17) 

Proof. First, for (<r, r) in S, 

il) r (o~,T) = I f(v)g(v + ra,v + rr)dv 
Je 

so that 

a 2k r 
M, Vv(^) = / f(v)R^g(v)dv , 

*>=o 7s 

where i?o-, r is the bidifferential operator given by 

l<h,i 2 ,...,i 2k <2d 11 12 l2k 

Hence 

d 2fc / d 2k 

Res(Cp(f®g),(3°) = (-) ^JLffiW) = V((^) |r= / 

= X /3 o( / f{v)Rfpg{v)dv) = [ f{v)Dfpg(v)dv 
Je Je 

where -D^o is the bidifferential operator given by 

D f3°9(v)= Yl ati,i 2 ,...,i 2k {P°)- ~ ~ (v,v 

l<ii,i 2 ,...,i 2k <2d 

where 



dv^dviz ...dvi 



2k(P°)= / Ph ■ ■ ■ Pi 2k Wf"\ T \^ 2 W ~ T f'M^r) , 



with the same convention as above. 

Of course these integrals are to be understood in the sense of distribu- 
tions, obtained by meromorphic continuation. □ 



12 



When k = 0, there is only one term, so D^o = c^o Id and it is possible 
to evaluate the constant c^o . 

Proposition 2.1. Assume that /3® + 0% + f3® = —2d, and assume that 0® £ 
-d-2N, j = 1,2 or 3. T/ien 



i? es (^(/),/3 ) = c (/3 ) / /(x)di 



where 



co(/3°) 



vr d r(2d) r( 



(2 V2)3^ r(f ) T{ z^) r{=%±) r(=f 



(18) 



(19) 



Proof. Let c be the stereographic projection, defined from E into the sphere 
5 of radius 1 in R d+l by 



2xi 



c(x) 



1+ b ■ 



i 2x d i 

\l+M 2 / 



The map c is conformal, and more precisely, for any tangent vector £ E 

2 



\Dc(x)£\ 



1 + Lr 



which implies 



/ f(a)da = 2 d [ f{c{x))(l + \x\ 2 r d dx 
Js Je 



for / any integrable function on S. Moreover, for any x, y in E, 

2\x - y\ 



\c(x) - c(y)\ 



(l + |af ) 1/2 (l + li/| 2 ) 1/2 



By the same change of variables, the trilinear form Cg is related to the 
trilinear form IC a on S studied in [3] through the relation 



K a (l ® 1 ® 1) = 2 Ql+a2+Q3+3 ^(/ ) 



(20) 



where ay = (3j + p for j = 1, 2, 3, and /o is the function on E x E x E given 
by 

/o0zi,X 2 ,Z3) = (1 + kll 2 ) ^(1 + |X 2 | 2 ) ^(1 + |X 3 | 2 ) — . 
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The left handside of ([20]) has been computed in [7j, and by a different method 
in [1] (see also [3]), so that 

r (f\ = (^\ 3d r(/gl + & + & + 2d)T(^± d )T(^± d )T(^± d ) 

'2^2 r ( fe+fe+^ r ( /Wi+rf ) r ( /3i+ff2+^ ^ J 

Hence, 

-V5Fn« r(^)r(^)r(^) 



Now 

/ / (x,x J x)da:= f {1 + \x\ 2 )-^ 1+Ms Ux = [ (1 + \x\ 2 y 

JE JE JE 



2d 



dx 



poo 

voliS^ 1 ) / (l+r 2 )- 2d r d - 1 dr = ir2 
Jo 



r( 



2 ■ 



T(2d) 

from which the lemma follows. □ 

In the general case (when k > 1), the integrals seem difficult to compute. 
Instead, in order to evaluate more explicitly the residue, it is possible to use 
the stronger invariance properties of the trilinear forms, which translate into 
covariance properties of the bidifferential operator -Dgo. 



3 Conformally covariant bidifferential operators 

Introduce the group of conformal transformations of E. A local transforma- 
tion <I> of E is said to be conformal if, at any point x where $ is defined, 
and for any tangent vector £, 

\D$(x)£\ = k(x)|£| . 

where k is a smooth strictly positive function, called the conformal factor 
of <E>. Classically, to any element of the group G = SOo(l,d + 1), one can 
attach a rational conformal action on E. If d > 3, then this group exhausts 
the group of positive local conformal diffeomorphisms (Liouville's theorem). 
The group G operates globally on the sphere 5 = S d of dimension d (see 
[21] ) and this action can be transferred to a (not everywhere defined) action 
on E by using the stereographic projection. It can also be realized as the 
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group generated by the translations, the rotations, the dilations and the 
symmetry-inversion i 

x 

\x\ 2 

To any element g € G, let n(g, x) be its conformal factor. Then k satisfies a 
co cycle relation, namely 

k(9i92,x) = K(g 1 ,g 2 {x))K(g 2 ,x) . 

A family of representations is associated to that cocycle. For A in C, define 

7r x (g)f(x) = K(g-\xy +x f(g- 1 (x)) , 

where we set p = i . These representations are the noncompact realization 
of the principal spherical series of SOq(1, d + 1) (cf [2Tj). 

Let A = (Ai,A2,As) € C 3 . A continuous trilinear form C on C^°(E) x 
C^°(E)xC^°(E) is said to be conformally invariant with respect to (ir\ 1 , ir\ 2 , 7r^ 3 ) 
if, for three functions fi,f2,f3 in C%°(E), for any g in a (sufficently small) 
neighbourhood of the neutral element in G, 

£(tt Ai ((7)/i ® TT X2 (g)f 2 ® n Xs (g)h) = ^(/i <8> /a ® /a) • (22) 

Recall the main result of [3], which was stated for the action of the 
group G on the sphere, but the situations are essentially equivalent through 
a stereographic projection. 

Proposition 3.1. Let (3 £ C 3 and assume that none of the conditions ([5]), 
([6]) is satisfied. Let A = (Ai,A2,Aa) be the unique element o/C 3 defined by 
the equations 

Pi = ~ Ai + A 2 + A 3 - p 

(3 2 = Ai - A 2 + A 3 - p (23) 
ft = Ai + A 2 - A 3 - p . 

i) The trilinear form Cp is conformally invariant with respect to tt\ x ® tt\ 2 

ii) Any continuous trilinear form which is conformally invariant with respect 
to ir\ 1 (g> 7T\ 2 71-;^ is proportional to Lp. 

By analytic continuation, the invariance is also valid for the residue of 
Lp at some pole (3°. At poles of second type, this invariance property can 
in turn be translated in a covariance property for the bidifferential operator 
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Definition 3.1. Let D be a bidifferential operator from C^ 3 (E x E) into 
C£°(E). Let A, /j,, v be three complex numbers. Then D is said to be covariant 
with respect to (tt\ ® tt^,tt u ) if for any functions f G C^°(E x E) and g in 
a (small enough) neighbourhood of the neutral element in G, 

Recall the following duality result for the representations tt\. 
Proposition 3.2. Let A G C. Then, for any functions <p and ip in C^ C (E) 

/ 7r\(g)(p(x)iJ)(x)dx = / y{x)it-\{g~ l )il){x)dx . 

J E JE 

This duality result links together (singular) conformally invariant trilin- 
ear forms and covariant bidifferential operators. 

Proposition 3.3. Let D be a bidifferential operator from C£°(E x E) into 
C^°(E). Let C be the continuous trilinear form defined for f G C^°(E) and 
g G C?{E x E) by 

C(f®g) = [ f(x)Dg(x)dx . 
JE 

Let A, fi, v be three complex numbers. Then the form C is invariant with 
respect to (irx,iT^,Tr u ) if and only if D is covariant with respect to (tt^ ® 
it v , tt-a)- 

Corollary 3.1. Let (3 = {/3\, 02, /3s) G %k f or some k G N, and such that 
none of the conditions ([5]) is satisfied. Let Dp be the associated bidifferential 
operator given by Theorem \2.4\ Let A = (Ai,A2,A3) be given by equations 
(f2Uj) . Then the bidifferential operator Dp is covariant with respect to (tt\ 2 ® 

This is a consequence of the expression of the residue Theorem 12.21 
together with the duality result (Proposition 13.2] ). The fact that none of 
the conditions ([5]) is satisfied amounts to the conditions 

A 2 ,A 3 ^-A: + N, A 2 + A 3 i -p-k-N . 

Covariant differential operators for the conformal group have been stud- 
ied intensively, and the following result was obtained sometimes ago by V. 
Ovsienko and P. Redou (see [IZ])- Recall the Pochhammer's symbol, for a a 
complex number and m G N 

(a) m = a(a + l)(a + 2) . . . (a + m — 1) . 
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Proposition 3.4. Let k be a nonnegative integer, and let A, p be complex 
numbers. 

i) Assume that A, p ^ {0, —1, —2, . . . , — {k— 1)}. T/ien i/iere exists a bidif- 
ferential operator which is covariant with respect to (vta®^, 7r\ + ^ +p+ 2k ) • 

ii) Assume moreover that X, p ^ {— p, —p — 1, ... , — p — {k — 1)}. Then 
the operator is unique up to a constant. 

(k) 

The operator D y x ' is explicitly described. The three fundamental bidif- 
ferential operators are A y , A z and the operator R defined for / in C°° (E x E) 
by 

pi d Vj dz j 

Then 

r,s,t,r+s+t=k 

where the c rs t are explicitly determined coefficients, depending on A, p and 
k, namely 

_(-iy- r fr + s + t\ (s + 1 

Crst 



2 r r\ V 1 ) (A + l) r 

r (24) 

r\t\ (X + p + r-s+ p) t - p (p + p + s + 2t) r _ p v ' 
when r <t, and for r >t, c rs t{\, p) = ct sr (p, A). 

Theorem 3.1. Ze£ /3 € /or some A; in N. Assume moreover that none 
of the conditions ([5]) are satisfied. Let X = (Ai,A2,As) the unique solution 
of the system J23J. T/ien (A 2 , A 3 ) i {0, -1, -2, (k - 1)} and 

^ = cpD^ (25) 

/or some constant cp. 

Proof. The conditions on 3 imply that A2, A3 4 — k + N, so that the con- 

(k) 

ditions for the defintion of D; are satisfied. Assume for a while that 3 
is such that A2, A3 ^ {0, — p, —p — 1, ... , —p — (k — 1)}. By the uniqueness 

(k) 

statement, there exists a constant cp such that Dp = cpD^\ 3 - By analytic 
continuation in the plane T~Lk, this equality remains valid on the domains 
were both sides are defined. □ 
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For instance, if k = 1, 

^ = - XTT A, + 2i?-- TT A z . (26) 

Theorem 13 , 1 1 however does not determinate the constant ca, and it seems 

quite difficult to test the operator D\^\ against the function /o as we did 
for the determination of the residue at a pole in %q. 



4 Bernstein-Sato identity 

A Bernstein- Sato identity (on the first parameter) is an identity of the form 

where (3 + 2 1 = (f3 1 + 2,/3 2 ,/3 3 ), B = B((x, y,z), 8 X , 8 y , 8 Z , (3) is a differ- 
ential operator with polynomial coefficients on E x E x E and depending 
polynomially in f3, and b is a polynomial in three complex variables. Such 
identities are known to exist (see [19], [20]), but are in general very difficult 
to find. It turns out that, in the case at hand, it is possible to find such 
identities. The proof uses in a crucial way the covariance property of the 
kernel In with respect to the conformal action of G on E. 

During the proof of some results, we will need to use the Euclidean 
Fourier transform. So it requires to extend the trilinear form to the Schwartz 
space S{E x E x E). This is merely routine. For the definition of invariance, 
one should formulate the condition in terms of the infinitesimal action of the 
conformal group. It is a classical computation (see e.g. [E]) and the action 
of the Lie algebra g = so(l, d + 1) or of the universal enveloping algebra is 
by differential operators with polynomial coefficients. Hence they operate 
on the Schwartz space S{E). Moreover the meromorphic continuation of la 
yields tempered distributions (cf jlO] for similar examples). Details are left 
to the reader. 

Lemma 4.1. Let M be the operator on S{E x E) given by 

Mip(y,z) = \y- z\ 2 ip(y,z) . 
Let A, /i be two complex parameters. Then M is an intertwining operator for 

Proof. Let g £ G, and ip E S{E x E), and assume that g is defined on a 
neighbourhood of Supp(ip). Then 

[Mo(-K X (g)m»(g))My,z) = \y-z\ 2 K{g-\yY +X K{g-\zY + ^{g- l {y),g- l {z)) 
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= \9-\y)-g- 1 {z)\ 2 K{g-\yY +x - 1 K{g-\zy^- 1 V {g- 1 {y\g- x iz)) 
= [n\-i(g) 7r^i{g)](Mip) (y,z) . 

□ 

Introduce now the Knapp-Stein intertwining operator. For v a complex 
parameter, let I u be the operator on S(E) given by 

W)(x)= [ \x-y\- d+ »f{y)dy . 
Je 

For 5if > 0, the integral is convergent and defines a continuous operator 
on S(E). It can be meromorphically continued to C, with simple poles at 
v = —2k, k € N. It satisfies the following intertwining property 

hv ° 7iV (5) = n-v(g) ° hu ■ (27) 

Now, for A, [i two complex parameters, form the operator 

N X:IX = I-2X-2 ® 1-2/1-2 M o I 2fl <g> I 2fi 

S (E x E) S{E x E) ^ S{E x E) !^^!z^ S (E x E) . 

For generic values of the parameters (A, fi), N\ tlJj is a well defined operator 
on S(E x E), which, by construction intertwines the representation tt\ ® ir^ 
and 7T A+ i g> tt^+i. 

Let A = (Ai, A2, A3) be a generic triple in C 3 , let (3 = (/3i, P2, P3) be the 
triplet associated to A through (|23|) . Observe that f3 + 2\ is associated to 
the triple (Ai,A 2 + 1,A 3 + 1). 

Consider the continuous trilinear form C on C^°(E) x C^°(E) x C^°(E) 
given by 

/2, / 3 ) = £{3+2! (fl ® iVA 2 ,A 3 (/2 ® / 3 )) • 
From the intertwining property of N\ 2 \ 3 

^(ttax (g)h,n\ 2 (g)f2, 7ta 3 (flO/3) = £/3+2! (ttaj (g)fi®Nx 2 ,\ a [tta 2 {g)f2®-K\ z {g)h}) 

= £/3+2i(7rAi(0)/l [tta 2+ i(5) ® 7r A 3 + l(ff)] O Nx 2 ,\ 3 [f2 ® /a]) 
= £ /3+2l (/ 1 0iVA 2 ,A3(/2®/3)) 

= £(h,f2, h) , 
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so that the form C is invariant w.r.t. (ir\ 1 ,ir\ 2 ,ir\ 3 ). By the generic unique- 
ness result on the invariant trilinear form (see [3]), the form C has to be 
proportional to Cp. Hence there exists a constant e = e(/3) such that 

(N X2M ) t (l f3+2l )=e((3)l f3 . (28) 

As we will see now, the operator (hence also its transpose) is a differ- 
ential operator on E x E, so that (|28p is indeed a Bernstein-Sato identity. 



Proposition 4.1. For X, fi £ C 2 , let E\^ be the differential operator on 
E x E defined by 

Ex,» = \y - z\ 2 A y A z 

d . ..A, , d 



3=1 3 3=1 3 

+ 20(20 + 2- d)A y + 2A(2A + 2- d))A z 



(29) 



8A V 3 9 

Its transpose F\ „ = E\ is given by 



F\» =\y~ z\ 2 A y A z 

d q d q 

+ 4(0 + 1) J2( z 3 ~ V3)q^ A v + 4 ( A + !) !>* - ^'fc A ^ 
+ 4(0 + 1)(0 + p)A y + 4(A + 1)(A + p)A z 



(30) 



- 8 (A + i) (P+ i)y:^ 

•7=1 J J 

The operator N\ tfl (for generic (A, fj,)) is a differential operator on E x E, 
and is given by 

N\,n = c(A, n)F x>tl , 



where 



= tt^ r(A)r(-A - i)r( M )r(-0 - 1) 



16 r(p - X)T(p + A + l)r(p - n)T(p + + 1) 
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Proof. Introduce the Fourier transform on E, defined by 

?m) = ho = [ e- i ^f(x)dx , 

JExE 

and extend it by duality to S'(E). The Fourier transform on E x E is defined 
accordingly. Observe that l v is a convolution operator with a tempered 
distribution, so that the Fourier transform of I v f is given by the product of 
the Fourier transform, i.e. 

n^m) = <*y)\z\- v fa) 

where 



(see e.g. [TO]). 

Next, as M acts by multiplication by a polynomial, the Fourier transform 
of Mip is given by 

T(M<p)(£, V ) = (-A f + 2R- 7]), 

where A is the Laplacian on E and R is the differential operator on S{E x E) 
defined by 

To prove the formula, it is enough to prove it for functions (p = f <8>fif, where 
f,geS(E). Now 



3 J \S)i 



so that 



A € (ir 2A /(0) = ier 2A A/(o 

d 



.7=1 °V 



+ 2A(2A + 2-d)|er 2A " 2 /(0 
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and 

9r u 

j=l 'J 
3=1 ^' 

+ 4A^|er 2A - 2 |r ? |- 2 ^ 2 (^^r/,)/(09W- 

3=1 

Let apart the factor c(2A)c(2/x)c(— 2A — 2)c(— 2/x — 2), the Fourier transform 
of N x>fJt (f (g> g) is given by 

- \e^f(0\v\ 2 9(v) +4A^e#(e)|r?| 2 5(r?) 

3=1 ^ 

- 2A(2A + 2 - d)f(£)\v\ 2 g(v) + 2Y,W 2 7ir^\ 2 p-(v) 

- 4A^3/(Oh| 2 |^(r ? ) " |e| 2 ^(e)^(r/) 

3=1 3=1 ^ 

3=1 

+ 4^]T |e| 2 /(6^# W - 2/1(2/1 + 2 - d)|£| 2 /(6^) • 

3=1 



Now use the classical formulae 

H«> fe7)«)=^f(«) 

a/({) = -If l 2 /(0 IviVGOK) = -a/K) 

to obtain the following expression for N\ jfJi (f<g>g) (up to the factor c(2A)c(2/u)c(— 2A- 
2)c(-2/i-2)) : 
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A y (\y\ 2 f)A z g + AXj^ JL(y jf )A z g 

3=1 Vj 

d 

+ 2A(2A + 2 - d)fA z g - 2 A y ( yj f)A z (z jg ) 

3=1 

~ 4A E ^^3) - ^A y (y 3 f)^-g 
3=1 Vj 3=1 3 

+ 4^ A y /— + 2^(2,u + 2 - d)A y / 5 

3=1 3 



The final expression for N\ ifl and (Nx^Y follows easily. □ 

As announced, E\ 2j \ 3 (being proportional to iV£ 2 ^ )is a candidate for 
a Bernstein-Sato identity for the kernel ip (the A's being related to (3 by 
(|23p). By brute force computation, the following identity is obtained. 



Theorem 4.1 (Bernstein-Sato identity). For (3 = (Pi, /3 2 , Pz) in C 3 , let Bp 

be the following differential operator on E x E 



Bp =\y - z\ 2 A y A z 



+ 2(/3 3 + A + d) - %)^-A 2 + 2(/3 2 + Pi + d) J>,- - z^—Ay 

3 = 1 Uj 3=1 3 

+ (p 3 + Pi + d)(p 3 +pi + 2)A Z + (p 2 + pi + d)(p 2 +Pi + 2)A y 
-2(/3 3 + /3 1 + d)(/3 2 + /3 1 + d)]T— — - . 

j=i °y3 oz 3 

Then 

Bplp +2l =HP)lp (31) 

where 

b(0) = (Pi + d)(Pi + 2)(/3i + p 2 + P 3 + 2d) (A + p 2 + p 3 + d + 2) . 
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5 Applications of the Bernstein-Sato identity 

The first application of the Bernstein-Sato identity is the computation of 
the residues of the distribution Cp along the plane %k by induction over k. 

Proposition 5.1. Let (3° be such that 0° + 0° + 0° = -2d -2k -2 for some 
k G N. Assume that £ -d - 2N (j = 1,2,3), and 0? / -2. Then, for 
f G S{E) and g G S(E x E), 

Re s{ C„^)Um) = (2k - - ^ - ff - JMO. 3° + 2,)(/^ ,) 

(32) 

Proof. For generic values of /3, from the Bernstein-Sato identitiy, 

= {tp,f®9) = ^{Bptp+z^ftog) = ^-{i p+2l ,f®Bpg) , 

and compute the residue at /3° on both sides. □ 

Let C/3 = Bp. Except for the change of parameters, this is nothing but 
the operator F\ jfl . 

Proposition 5.2. 



C P = B%= \y-z\ 2 A y A 



d 



+ 2(a + p 2 + d + 2) - yj)*- A 



dz~ v 

d 



+ 2(0i + ft + d + 2) Yftj - Zj ) A Aa 
+ (A + 02 + 2d) (0i + /3 2 + d + 2)A y 
- 2(0! + /3 2 + d + 2) (01 + /3 3 + d + 2) ^ — — - 

j =l a yj az j 

+ (0! + 03 + 2d) (0i + 03 + d + 2)A, . 

To write an expression for the residue at a pole in T-Lj., where k G N, let 
use the following convention: for (3 = (01,02,03) arid G N, let 

/3-(2fc)i = (0i-2fc,0 2 ,0 3 ) 

Now, for (3 G "Ho, define the differential operator Eg on E x E by 

Cp 0) = U, Cp k) =C H o...oC m 
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Theorem 5.1. Let (3° G H , and let k £ ~N. A ssume that j3j ^ —d — 2N 
(j = 1, 2, 3) and $ (£ {0, 2, . . . , 2k - 2}. T/ien 



Res(C p ,0> - (2&)i)(/® 5 ) = c fc (/3°) / f(x)(C$g)(x,x)dx (34) 

Cfc ^°) = liiin — * ^ Co(/30) • 

16**1 (p)* ( _^ )jk( _^_ p + 1)jb 

The Bernstein-Sato operator can also be used to describe a family of 

(k) 

covariant bidifferential operators. Let A, \i be in C. For k G N, let be 
the bidifferential operator defined by 

F ^f( x ) = Fx+k-i^+k-i o ... o F x ^f(x,x), (35) 



where F\^ is the differential operator on E x E defined by (|30lh 

Theorem 5.2. Lei A, [i in C, and A; € N. Then the operator F^ 1 is confor- 
mally covariant with respect to (ir\ <E) ^^,^x+^+p+2k)- 

Proof. Recall that the operator F\ :fl is covariant w.r.t. (7r A (g>7r M , n\+i®TTn+i)- 
So, by induction, F A+fe _i ifl+fe _io- • -oF x ^ is covariant w.r.t. (7r,\®7r M , 7r A+fc (g> 
TT^+k)- Now the map 

where /(x) = f(x,x), is covariant w.r.t. (7r A+fc ® Ti^+fc, 7T A+(U+p+2 fc)- The 
assertion follows. □ 

For k = 1, one gets 

F fi = + !)(/* + P) A 2/ " 8 ( A + 1)(A* + 1)^ + 4 ( A + + P)A, 
to be compared with (1261) . 



6 Final remarks 



1. The construction of the covariant differential operator iV A „ admits a 
natural generalization. Let r be the standard representation of G and r' its 
dual representation. Choose highest weight vectors v and eft (with respect to 
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some suitable ordering) for r and r', respectively. Then up to some constant, 
the multiplication by the matrix coefficient 

9X,92 ' — > (-r(gi)v,T'(g 2 )4>) 

coincides with M from lemma H~T1 Upon replacing r by another irreducible 
finite-dimensional representation of G one obtains a multiplication operator 
that intertwines the tensor products of general (i.e. not necessarily spherical) 
principle series representations (see [1]). Using appropriate intertwining 
operators one obtains an associated differential operator. We would also like 
to point out that the operators used by Oksak in his work [16] on invariant 
trilinear forms for G = Sl 2 (C) — Spin(3, 1) are of this type. See also |12j . 

2. A byproduct of Proposition 12.11 is that the residue at a point (3 G Hq 
vanishes identically if ~^\~ d G — N, i.e. if f3\ G — d + 2m, m G N. For such a 
value, observe that 

x P2 + P3 , -Pi -2d , 

Ai = — - — + p = +p=-m . 

Now for m > 1, there exists a nontrivial differential operator i? m (see [2]) 
on C^?°(-E') which is covariant w.r.t. (7r_ m ,7r m ). Now let 

3 = (01 - 2m, /3 2 + 2m, /3 3 + 2m) , 

so that 

A = (m,A 2 ,A 3 ) . 

As j3\ + /?2 + /?3 = —2d + 2m, /3 is no longer a pole. So the form is well 
defined, and the form 

{fl, f2, h) 1 > £p{Rmfl, f2, /3) 

is invariant with respect to (7r_ m , ttx 2 ,7T\ 3 ). However, the relation Xi + X 2 + 
A3 = — p guarantees that the form 

(/1,/2,/s)'— > / fx(x)f 2 (x)h(x)dx 
Je 

is invariant under (-7r^ 1 , 7r^ 2 , -k\, a ) . So, for Ai = — m and (A 2 ,A3) generic, 
we have produced two (linearly independant) trilinear invariant forms on 
C^°(E)xC^°(E)xC^(E) w.r.t. (vr_ m , vr A2 , 7r A3 ) . Although we won't develop 
these aspects here, the same remark can be used to produce, for specific val- 
ues of A, two (linearly independant) bidifferential operators covariant under 
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the same actions of the conformal group. Notice that this is in concordance 
with the results and the philosophy of [13] and |14j . 

3. The relation between our formula for covariant bidifferential operators 
(|35p and the formulae obtained in [17] or in [14] is still to be investigated, 
and the coefficients which relate them seem to be important. In the classical 
setting (i.e. for the original Rankin-Cohen operators acting on the upper 
half-plane), much effort has been devoted to understand the structure of this 
family of operators (see [23], [5], [22], [6], [IB]). We hope that our realization 
of these operators will add to the understanding of the family of generalized 
Rankin Cohen operators. 
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